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Abstract 

For a matrix *-algebra B C C™^™, denote by ^ := Sym"(B) the matrix *-algebra 
consisting of the elements in the n-fold tensor product B'^" that are invariant under per- 
muting the n factors in the tensor product. Examples of such algebras in coding theory 
include the Bose-Mesner algebra and Terwilliger algebra of the (non)binary Hamming 
cube and algebras arising in SDP-hierarchies for codes using moment matrices. We give 
a computationally efficient block diagonalization of A in terms of a given block diago- 
nalization of B. As a tool we use some basic facts from the representation theory of the 
symmetric group. 

Keywords: block diagonalization, semidefinite programming, Terwilliger algebra, as- 
sociation scheme, representation theory, symmetric group. 

1 Introduction 

A matrix *-algebra is a set A C c^x" of matrices that is closed under addition, scalar 
multiplication, matrix multiplication and A ^ A* (taking the conjugate transpose). It is a 
classical result that any such algebra can be brought into block diagonal form. That is, there 
exists an isomorphism 

t 

0:>1^0CP«^P' (1) 
1=1 

of *-algebras, meaning that (/> is a linear bijection, (j){AB) = (/>(^)(/)(i?) and 4>{A*) = </>(^)* 
for all A,B e A. 

In this paper we will be concerned with constructing in a computationally efficient way 
such block diagonalizations (p. The motivation comes from semidefinite programming, where 
block diagonalization has proven to be a valuable tool since it can be used to reduce the 
complexity of semidefinite programs having a large group of symmetries. 

Semidefinite programming is an extension of linear programming that is both very general 
and at the same time can be performed efficiently (up to given precision), both in theory and 
in practice. A (complex) semidefinite program is an optimization problem of the form 

maximize {C, X) subject to X ^ and {Ai, X) = hi for i = 1, . . . ,m, (2) 
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where X £ c^x" ig Hermitian matrix variable, C,Ai,..., Am G C"^"" are given Hermitian 
matrices and bi,...,bm are given real numbers. Here X ^ means that X is positive 
semidefinite and {C,X) := trace X*C is the trace product. Linear programming can be 
viewed as the special case where the given matrices are diagonal. In the literature, it is more 
common to consider semidefinite programs where the matrices are real valued and symmetric 
rather than Hermitian. However, complex semidefinite programming is easily reduced to the 

real case by mapping each Hermitian matrix A to ReA^) ' ^ ' 

In recent years, many results have been obtained using semidefinite programming, where 
the underlying problem exhibits a large group of symmetries (see [7] for an overview) . These 
symmetries can be exploited to significantly reduce the computational complexity of the 
semidefinite program at hand. Indeed, let ^ be a matrix *-algebra containing the given 
matrices C,Ai,..., Am- For example, A may be the set of matrices invariant under a group 
of symmetries of the underlying combinatorial problem, or the matrices C, Ai, . . . , Am may 
belong to the algebra A associated with an association scheme (or coherent configuration). 
Then the matrix-variable X can be restricted to A without changing the optimum. When an 
explicit *-isomorhism (j) : A ^ ^^C^*^P» is known, the semidefinite program can be reduced 
to a semidefinite program in terms of the smaller matrices from C^'^P' using the fact that 
a *-isomorphism preserves positive semidefiniteness. When the algebra A is commutative, 
for instance when A is an association scheme, the optimization problem is hence reduced to 
a linear program in a number of variables equal to the dimension of A. 

When applying block diagonalization to semidefinite programs, it is necessary to have a 
*-isomorphism available that can be effectively computed. When A consists of the matrices 
invariant under the action of a group, the theory of group representations is the tool for 
deriving an explicit block diagonalization of the invariant algebra A, see [7]. 

In [U m [H] semidefinite programming is used to obtain bounds for error-correcting codes. 
There the symmetry reduction is essential, reducing the size of the matrices from exponential 
to polynomial in the size of the input. The underlying combinatorial problem involves strings 
of length n over a finite set of symbols and the main symmetry comes from permuting the n 
positions in a string. This leads to the consideration of algebras of a specific form. 

Given a matrix *-algebra B and an integer n, the algebra B®^ admits an action of the 
symmetric group 5„ by defining a{Bi^- ■ -(^Bn) := i(i)'^' ' ■'5?'-Bo--i(n) for . . . , Bn € B, 
a (z Sn and linearly extending this action to the whole of B®^. The set of symmetric tensors 
Sym"(i3) := {A £ B®"" \ aA = A for all a G Sn} is a matrix *-subalgebra of 

In this paper we derive an explicit block diagonalization of A := Sym"'(,B) in terms of a 
given block diagonalization for B. The resulting block diagonalization can be computed in 
polynomial time from the block diagonalization of B. When B and its block diagonalization 
are defined over the reals (rationals), also the obtained block diagonalization of A is defined 
over the reals (rationals), up to a scaling by square roots of rationals in the latter case. 

In Section ??, we will consider some examples, including the Terwilliger algebra of the 
binary- and nonbinary Hamming cube taking respectively B = C^^^ and ;B = C^^^©C^^^. 
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2 Preliminaries 



2.1 Notation 

We denote the natural numbers by N := {0, 1, 2, . . .} and the positive integers by Z+. For 
k € Z+, we denote [k] := {1, . . . ,k}. For finite sets S,I and a word a € , the weight of 
a, denoted w{a) G N'^ is given by w{a)s := |{i G / | Oj = s}| and counts the number of 
occurences of each element s € S in the word a. 

Let V he a vector space (say over the complex numbers) with basis vi, . . . ,Vk- A basis 
of the n-fold tensor product 1^®" is given by (fa)ae[fe]"; where f a = ^^ai ® ■ ■ ■ (X" Va„- The 
subspace Sym"(y) C y'^'^ consisting of the symmetric tensors has a basis indexed by all 
decompositions ^ G of n into k parts: := ^a£[k]",w{a)=^l'^a.■ 

2.2 Representation theory 

In this section we recall some basic facts from the representation theory of finite groups. 
These results can be found in most textbooks on group representations such as the book by 
Fulton and Harris [2]. Given some finite set S and a finite group G acting on S, we derive a 
*-algebra isomorphism between the set of complex S x S matrices that are invariant under 
the action of G, and a direct sum of full matrix algebras. This isomorphism is expressed in 
terms of representations of G. 

Let y be a vector space. All vector spaces that we consider will be finite dimensional 
over the field of complex numbers. By Gh(V) we denote the set of all invertible linear 
transformations of V to itself. Let G be a finite group. Then V is called a G-module if there 
is a group homomorphism p : G ^ GIj{V). That is, there is an action of G on 1/ such that 
g{c ■ V + c' ■ w) = c - g{v) + c' • g{w) for all g G G, c,c' € C and v,w G V. 

A G -homomorphism from V to W, is a linear map (p : V ^ W that respects the action 
of G: 4>{g{v)) = g{(j){v)) for all v G V and g G G. If <j) is a bijection, then V and W are 
isomorphic (as G-modules) and we write V = W. The set of homomorphisms from y to 
is denoted by Home {V, W) and we denote EndG(^) := Home {V, V). 

Let (•,•) be a G-invariant inner product on V, that is (x,y) = {g{x), g{y)) for all 5 € G 
and x,y GV. Such an inner product exists: take {x,y)Q := X^^gc idi^)^ div)) inner 
product (•, ■). With respect to this inner product, the algebra EndG'(^) becomes a *-algebra, 
where A* is the adjoint of ^ € EndciV). 

A module V is called irreducible if V has exactly two submodules: {0} and V itself. If 

C y is a submodule, then also W-^ := {x £ V \ (x, y) = for all y E W} is a submodule. 
Hence it follows that V can be decomposed into pairwise orthogonal irreducible submodules. 

Theorem 1. (Maschke's Theorem) Let V he a G-module. Then V has an orthogonal decom- 
position 

k mx 

^ = 0^A, ^A = 0n^ (3) 

A=l i=l 

where = for all i, j = 1, . . . , mx and the Vjl^ (X = 1, . . . , k) are pairwise nonisomorphic 
irreducible submodules ofV. 
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The set of homomorphisms is controlled by Schur's Lemma. 

Theorem 2. (Schur's Lemma) Let V and W be irreducible G -modules. Then either there 
is an isomorphism (f) : V ^ W and Home {V, W) consists of the scalar multiples of <j), or 
Home {V, W) consists of the zero map only. 

Together, Schur's Lemma and Maschke's Theorem imply a block diagonalization ^xidciy) ? 
©A=i C™^^"^^. To make this explicit, let f/i, . . . , [7^ be a complete set of irreducible submod- 
ules of V, say Ux = for A = 1, . . . , k. Let ex^Ux\ {0}, Wx ■= {Aex \ A € EndG(l^)}, and 
let Bx be an orthonormal base of Wx. Then we obtain a block diagonalization of EndG(y) 
as follows. 



Theorem 3. The map 

V':EndG(F) ^ 0C™^x'™A^ (4) 



k 



A=l 

k 



A ^ ®{{Ab,b') 



'h,b'eBx 
A=l 

is a *-algebra isomorphism. 

Proof. For each A, the map A i— > {{Ab, b'))^ b'eBx ^ *-algebra homomorphism since it maps 
A to its restriction on Wx written as a matrix with respect to an orthonormal base of Wx- 
Therefore ip itself is a *-algebra homomorphism. 

To show injectivity of ip, suppose that 'ip{A) = 0. Consider an arbitrary component 
in Maschke's decomposition ([3]) and let B : Ux ^ he a G-isomorphism. Then B can be 
viewed as an element of EndG(y) by setting Bx := for all x G U^. Since ip{A) = it 
follows that ABex = and hence 

A-V^ = AB ■Ux = AB ■ CGex = CG ■ ABex = {0}. (5) 

Since this holds for all A and alH, A is the zero map. Surjectivity follows since 

dimEndG(l/) = dimHomc 00n^00^; (6) 



M i=i 



nix ™M 



= EEEdimHomG(n\^;i) 

i=l j=l 

= E"^A> 
A 

shows that both the dimension of FjnddV) and the dimension of ©^ C»"AxnT-A equal m|. 
Hence ip is indeed a bijection. □ 

The theorem shows that once a nonzero element ex in each irreducible submodule (up to 
equivalence) is identified, the block diagonalization of EndG(y) can be computed explicitly 
once we can evaluate the inner products {AA'ex,A"ex) for any ex and any A, A', A" from a 
suitable basis of EndG(V')- 
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3 The algebra Sym"(C^'^^) 

Consider the set [p]" of words of length n with symbols in the alphabet [p] = {1, . . . ,p}. For 
example, taking p = 2 gives the set of all binary words of length n written with symbols 1 
and 2. The symmetric group Sn acts on [p]" by {cra)i := a^-i(^i^ for all a G [p]", a e Sn 
and i = 1, . . . ,n. This indiiccs a linear action of Sn on C'^' by defining {(Jx)a := a;^— 1„ for 
X € C^^^ . Similarly, there is an induced linear action of Sn on the set of [p]" x [p]" matrices 
M, by defining (aM)„,(. := M^-i^^^-i^. 

In this section, we give a block diagonalization of the the matrix *-algebra 

A = {Ae C[P]"^W" \aA = A for all a € Sn} (7) 

of all S'„-invariant matrices in Ct^' ^^^^ . That A is indeed a matrix * algebra follows from 
the facts that Sn acts linearly on C^p^^Ip]" , a A* = {a A)* and a(AB) = {aA){aB) for all 
a eSn and A, B e C^^" ""^^^ . 

For any two words a,b e [p]", define D{a, b) G N^^^ by 

{D{a,b))ij := \{k \ ak = i,bk = (8) 

Clearly, for a, b, a' , b' G [p]" we have D{a, b) = D{a\ b') if and only if a' = 7r(a),b' = 7r(b) for 
some TT e Sn- Let 

P{n,p) := {D e NP^'P \ l^Dl = n}. (9) 
For D G P{n,p), let G v4 be given by 

(^I3ja,6 — < . (10) 

10 otherwise. 

Observe that Aj^ = Aj^t. 

Proposition 1. The matrices A^ with D G P{n,p), form a basis for A (as a complex linear 
space) and the dimension of A equals {^''p2_i^)- 

Proof. The matrices A^ are nonzero matrices with disjoint support and hence linearly inde- 
pendent. For any matrix A G A, the value of an entry A^^b only depends on D{a, b) and is 
therefore a linear combination of the matrices Ar,. □ 

Since A is closed under multiplication, there exist numbers j^, L, M, N G P{n,p) such 
that AlAm = YliN M-^i^ ■ Although we will not need the numbers c^}^, we mention the 
following fact for completeness. 

Proposition 2. The numbers are given by 

V 

B r,t=l 

where the sum runs over all B G f^P^P'^P that satisfy Y2k^r,s,k = -^r-.s? ^k^k,s,t = Ms,t O'l^'d 
J2k ^r,k,t = Nr,t for all r, s, t G [p]. 
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Proof. For given words a, c G [p]" with D{a,c) = N, the number equals the number 
of words b € [p]™ with D(a,b) = L and D{h,c) = M. Fixing tlic number -Br-,s,t of positions 
i = 1, . . . ,n for which a-j = r, hi = s,Ci = t, the number of feasible b equals unless Br,s.k = 
^r,s-i X^fe -Br,fc,t = -/Vr,t and -Bfc^s^t = Ms^f for all r, s, t. In the latter case the solutions b 
are obtained by partitioning each set Sr,t '■= {i \ = r^Ci = t) into subsets of size Br^s,t for 
s = 1, . . . p to choose the entries of h on Sr^t- D 

3.1 Block diagonalization of Sym"(C^'^P) 

The action of Sn on C'^' gives C^^l the structure of an S'n-module. Identifying matrices 
in C^^l ^[^1 with the corresponding linear maps in GL(Ct^l ), the subalgebra A is identified 
with EndG(C[P] ). Using the representations of the symmetric group, we will find an explicit 
block diagonalization of A. 

A partition A of n, written A h n, is a sequence of nonnegative integers Ai > A2 > • • • > An 
with n = Ai + A2 + • • • + A„. Partition A is said to have k parts if exactly k of the numbers Aj 
are nonzero. A Ferrers diagram of shape A is an array of n boxes, taking the first Aj boxes 
from the i-th row of an n x n matrix of boxes. The j-th box in row i is referred to as the box 
in position We also number the boxes from 1 to n according to the lexicographic order 

on their positions: the box in position has number Ai + • • • + Ai_i +j. The dual partition 
A* of n gives the column lengths of a Ferrers diagram of shape A: A* = |{j | Aj > 

A tableau t of shape A is a filling of the boxes in the Ferrers diagram with integers. 
Here all entries will be from the set [p\. The number in box k is denoted by t{k). Fixing a 
partition A h n, each element a of [p]" is identified with the tableau t given by t{k) := o^. 
For an element a E Sn and a tableau t, we define {at){k) := t{a^^{k)). A tableau t is called 
semistandard if the entries are increasing in each column and non-decreasing in each row. 
We denote the set of semistandard tableaux of shape A with entries from [p] by ^,p- 

Given a partition A l- ra, we can associate with A two subgroups of Sn- The group C\ 
consists of the permutations of boxes within the columns of the diagram, and the group R\ 
consists of the permutations that permute the boxes within the rows. Two A-tableaux t, t' 
are called (row-) equivalent, written t ^ t', when t' = irt for some vr G Rx. 

Recall that the complex vector space y = is an S'n-module. We denote hy x"' ^ ^ 
the standard basis vector corresponding to the word a. Then the action of Sn on V is given 

by 

for all a Sn and S C for all a S [p]". 

Let A be a partition of n into at most p parts. Define the tableau tx by filling the positions 
in the i-th. row with z's. Given any tableau t of shape A, we define et EV hy 

et:= E sgn(a)Ex'^^*')- (13) 
Define S^ := CSn ■ ca, where ex := et^ = Eo-eCx sgn((T)x'^^*^^. 
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Theorem 4. The S'^, where A runs over all partitions of n into at most p parts, form a 
complete set of pairwise non-isomorphic, irreducible submodules ofV. 

The irreducible modules are called Specht modules. The proof is standard and can be 
found for example in [5]. 

Let D G P{n,p) and let A = D^l. If A is nonincreasing, we view it as a partition of n and 
make a tableau t{D) of shape A and weight ^ = Dl as follows: the i-th row of t{D) contains 
Dj^i symbols j, and we make the rows of t{D) non-decreasing. Observe that D(t(D), tx) = D. 
We have the following lemma. 

Lemma 1. Let A = Ad & A and let X' := Then 

Aex = h^^ '^^' = ^ (14) 
I otherwise. 

Proof. Recall that we identify [p]" with the set of all tableau of shape A and entries from 
[p\. For any such tableau t, we have Ax^ = Ylis\D(s t)=DX^ • Since D{s,tY\ = w{t), this sum 
equals zero when w{t) ^ A'. In particular, 

4 _ if A / A', 

''^'[T.s^tmX' ifA = A'. ^^^^ 
The lemma now follows by writing out the definition of Aex: 

Aex = sgn(a)^x''^*^^ (16) 

= Yl sgn{a)a{Ax^^). 



□ 



The lemma implies that the linear space Aex is spanned by the vectors ct (where t S [p]" 
is a tableau of shape A). The following theorem selects a subset of the tableaux to obtain a 
basis. 

Theorem 5. The et, t € T^ p constitute a basis o/span{ef | t £ [p]" is a tableau of shape A}. 
The proof is standard and can be found for example in [5] . 

In general the basis {et} {t G Tx,p) of Aex, is not orthonormal. However, let B := {et)t be 
the matrix with the et as columns and let Gx '■= B be the Gram matrix of the ej. Take 
a Cholesky decomposition {Gx)~^ ='■ RxRx^ then the columns of BRx form an orthonormal 
base of Aex- 



Theorem 6. The map 



ij-.A ^ QC"^^™^ (17) 

A 

A ^ ®Rl{{Aes,et)),^^Rx 



is a *-isomorphism. 
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Proof. This follows directly from Theorem [3] since the et, t Tx^p form a basis of the space 
Aex, where the ex generate the irreducible submodules of V. □ 

Remark 1. Although the vectors et have length exponential inn, we can compute their inner 
products efficiently (see next section). This implies that we can find the Gram matrices Gx 
efficiently. However, in applications to semidefinite programming, we only need ^ to preserve 
positive semidefiniteness. Hence we can neglect the matrices Rx and obtain a bijection 

ij':A ^ 0C'"^x'"A (18) 

A 
A 

that preserves positive semidefiniteness. 



3.2 Computing ip 

In this section, we wil show how the map ip can be computed efficiently. That is, given a 
partition A of n into at most p parts, and given semistandard tableaux s, i G 7a,p of shape A, 
we compute for every D G P{n,p), the inner product {ADes,et). 

FovAe cIpI'^IpI" the linear map A ^ C given hy Ad i-^ {Ad, A) is conveniently expressed 
using polynomials in the p'^ variables Xij, i, j = 1, . . . ,p: 

fA:= ^""{^D^A), (19) 

DeP(ra,p) 

where the shorthand notation used. 

Define for A; = 0, 1, ... ,p the polynomials Qk by setting Qq{x) := 1 and 

(Xl^l ■■■ Xl^k^ 

Qk{x) ■= k\ det 



(20) 



\Xk,l ■ ■ ■ Xk,k^ 



for k = 1,2, ... ,p. The polynomials Qk have at most p\ terms and can, for fixed p, be com- 
puted in constant time. Given a partition A h n into at most p parts, define the polynomial 
Px by 

n 

i=l 

For fixed p, computing Px can be done in time 0(ndim(^)). 
Proposition 3. Let X\- n be a partition with at most p parts. Then 

Px = fe,el- (22) 
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Proof. Writing out the definition oi Ad and ex we obtain 

x"" ■ ^AD,exel) = E E sgn(p)sgn(T)(Ai, V,,,*, (23) 

D&P{n,p) D&P{n,p) p,TeCx 

= E ■ E Sgn(^)(^^)<^«A,*A 

l?eP(n,p) <t6Cx 

= |Ca| E sgn((7)x-°(^*^'*^). 

Here we used the substitution a := p so that 

sgn(p)sgn(T)(^£))pt;^^T-tA = sS'^{^){-^D)atx,tx- 

(24) 

Any a E Cx corresponds to an n-tuple (cti, . . . , cr„), where crj G S'^t is the restriction of a 
to the z-th column. Observe that {tx)ij = i and {atx)i,j = '^J^i'i) and hence 



\ * 



^D{atx,tx) 



(25) 



Using this, we obtain: 



E sgn(^)^''^'^*^'*^^ = |CAin(E sg^C'^.OlIs^^W.^^ ^^^^ 
a-eCx j=i o-jeSx* 2=1 

n 



A- 



□ 



Since we can perform multiplication in A in polynomial time, we now have a polynomial 
time algorithm to compute {Acf, Cfi) = (^Aj^^^^, ^_^-^A{D)Ad^^^^^^, exej^"^ for any t,t' G Tx^p and 

D G P{n,p), hence obtaining the map ^ explicitly. 

Below we show how we can speed up the computation, by avoiding the computationally 
costly multiplication in A. 

For i,j G [p], let Ai^j G >A be defined by 

I 1 if there is a A; such that = i,hu = i and ai = hi for sill ^ k, , , 
{Ai^^)ah ■■={ k , k J I I ^ , ^27) 

lO otherwise. 

For i,j G [p], let Eij be the [p] x [p] matrix with {Eij)ij = 1 and all other entries equal 
to 0. 
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Proposition 4. For any D € P{n,p) we have: 

AdA^, = iDk,j + l)AD-E,,,+E,j, (28) 

k\Ak,i>0 

Ai-.jAD = Yl (A,fc + l)vlD-£;,.,+£;,,,. 

k\Aj,k>0 

Proof. To prove he first line, consider two words a,b £ [p]". The entry (A/jAj^j)^^^ equals the 
number of c € [p]" such that D{a, c) = D and b is obtained from c by changing an i into a j in 
some position h. If = k, then D{a, b) = D — E/^^i + -Efej- When D(a, b) = D — E^^i + -Efcj-, 
the number of possible c equals -D(a, = Dk,j + 1. 

For the second line, observe that Ai^jA^ = (A^jAj^i)'^ . □ 

Given a decomposition fi = (^i, . . . ,/ip) of n (^j > 0), let := ^^j^;^ j+.^+^p^^.p- 
Proposition 5. If D G P(n,p) is lower triangular with /j, := D^l, then 

u^^A = (li n A'i'i)^^- (29) 

^i>j J j=l i=j+l 

Proof. For given i ^ j and D E P{n,p) with Dj j > 1 and Dj^fc = for all k ^ j, Proposition 
(HI) gives 

Ai^jAo = (Dij + 1)Ad^e,,,+e,,, ■ (30) 

Let (ii, ji), . . . , {im,jm) be m pairs of indices with ik 7^ jfc for all k and > ji whenever k < I, 
and let si, . . . , Sm be nonnegative integers. Let D := + • • • + fJ.pEp^p + YlT=i ^ii-^ikJk ~ 

^jkJk) suppose that D is nonnegative. Then from (f30l) we obtain 

^?:->..---A%.^M = --!---«i!^^' (31) 
by induction on m. This implies the statement in the proposition. □ 
Let us denote 

p d ^ d 

^i^i y^.^i'S-^ , dj^j ■= y^^s.j Q- • (32) 

s=l "" s=l ' 

In terms of polynomials, Proposition (j3]) now gives 
Proposition 6. For i ^ j and a matrix A, we have 

fAi_jA = di^jfA (33) 
fAA,^j = d*^jfA- 
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Proof. We have: 

D D k\Dk,^>0 

k\Dk,r>0 D' 

Here we used the substitution D' := D — Ek^i + -E^j. The proof for fAi^jA is similar. □ 



Theorem 7. Let t',t" be semistandard X-tableau and let D' := D{t' ,tx), D" := D{t'\t\). 
Then 

x^{Ane,,e,.) = q-J\{D'^JDl^r)-\ (35) 

DeP{n,p) i>j 

where 

Q=UIl o ^'a. (36) 

j=i i=j+i 

Proof. We want to compute 

{ADet',et") = X]^^ (^AD,An"exelAl,'j = fA^^E^A^,- (37) 

D D ° 

By Proposition ([5]) we obtain 



p— 1 p p— 1 p 

.-D" .Vt -r-r .D' 



\{{D[,^Dl^).Ao>:ExAl, = J] 11 ^^^•FAi5A/A(n 11 A-I'j)^ (38) 

i>j j=l i=j+l j=l i=j+l 

p—l p P^l P 

= n n -^Sain n •^Sj)'^ 



j = l i=jr- + l j = l i=j + l 



Using Proposition ^ we derive 



ll(Di,lDlJ) ■ fA,„E,Ar, = li fl m^,f^<d*^,f'^.*) o P,. (39) 



D' 

i>j j=l i=j+l 



□ 
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4 The general case 

Let Bhe a complex matrix algebra with basis Ri, . . . , Rg, and let ^ := Sym"(S). Recall that 
we have a basis {i?^ | G N^} given by i?^ := J2x&[s]"\w{x)=^^ ^'Li^^xr 
Assume that a block diagonalization 

t t 
<^ : ^ ^ c^xfS 5^0 d^iiB) (40) 
i=l 1=1 

is given. In the following, we will describe a block diagonalization of A in terms of (f) and 
with respect to the basis {Rn}- 

First observe that the isomorphism : B®'^ ^^t^^f^Pixpi'^^n gjygg isomorphism 
A Sym"' {®j^iCP*^P*) by restriction to A. Now observe, that by removing multiple copies 
of identical blocks in Sym"(©*^jC^»xp*), we have an isomorphism 

Sym-(e*=iC^''XfO ^ e^eN* ® Sym^^ {CP^^'p^) (41) 

Indeed, let {Ej^i^i^ \ k,l = 1, ...,m} denote the standard basis of C". This gives a basis 

{Ed \ D e n;^^"*} given by Ed :— '^x&{[m]x[m])'^,w{x)=D -^x- Similarly, the standard basis 

{E(^i^f,^i I i = l,...,t,k,l = l,...,pi} of ©l^iCf*^*'' gives a basis {Ed \ D = Di ® ■ ■ ■ ® 

A} 'of Sym'^(e*^iCf^xfO where Ed := Exe([pi]xbi]u-uMxM)"|«,(x)=Die-eA- Then the 
isomorphism is given by 

Ad ^ e^SiTDa,^, ■ ■ ■ Voa,/.* ®l=i ^A- (42) 
This yields an isomorphism 

t 

^■.A^ (g) Sym'^^ (CP» XK ) _ (43) 

To make this map explicit in terms of the given bases, we use the following lemma. 

Lemma 2. Let V he a vector space with bases {ui, . . . , Um} and {vi, . . . , Vm} that are related 
by Vi = Ylj Introducing indeterminates xi, . . . , Xm, the bases {u^)^^fim, and {vy)y^mm 

o/Sym"'(F) are related by 

Vv = ^U^- (xiCi,! H 1- XmCl,mY^ ' ' ' {xiCm,l H 1" XmCm,mT' [x"]- (44) 
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Proof. We compute J2iym'" ^'^^i- 



(m 
E 
i=i 

= E 



j;,z£{l,...,m,}" 

(m 

m \ / m \ 



Application of the lemma yields the following result. 
Theorem 8. The map ip is given by 

t 



□ 



u tj. D={Di,...,Dt) i=l 

where l^Dil = fXi. The coefficients yo are given by 

(t V] / s \ {^j)fe,i\ 

n n [Y.'^^{^^^mKl\ [^1- (47) 
i=ifc,«=i \j=i / / 

Proof. Expressing the basis {^{A,^) \ v G N^} in terms of the basis {^d} using Lemma [5] 
yields the claimed result. □ 

Using the block diagonalization for each of the algebras Sym'^' (C'^'^^*) as described in 
section 13. H we obtain a block diagonalization 

t 

^^0 (8)C^^»-»"^^»-» (48) 

M Ai,...,At i=l 

of A 



13 



5 Examples from coding theory 
Terwilliger algebra of the binary Hamming scheme 

The algebra A := Sym"(C^^'^) is referred to in the hterature as the Terwilhger algebra of the 
Hamming scheme. The matrices D G P(n,2) are usually indexed by 3 parameters i,j,t by 
setting 

'n + t-i-j i-t\ .( 



Here < f < and i + j — t < n. In [6] an explicit block diagonalization was given and 
used to compute semidefinite programming bounds for binary codes. Here we show that our 
method gives the same block diagonalization. 

Since p = 2, the partitions X = [n — k, k) are indexed by = 0, . . . , [^J . Give k, the 
semistandard tableaux of shape {n — k,k) are indexed hy i = k, . . . ,n — khy placing a 2 in i — k 
of the last n — 2k boxes of the first row and in all k boxes of the second row of the tableau. 
Lets denote this tableau by t^^i and denote et^- by ek,i- We have D{tk,i,tk,k) = {^Ikk)- 
Hence for give 

It is easy to see that (^2^1 (2^2,22^1,1 — 2:2,1X1^2) = ^1^2(^2, 22^1,1 — 2:2, 1X1^2) = 0. Hence we 
obtain: 

2^ 

{ADek,i,ekj) = (• _ _ fe), (2^2,2x1,1 - xi,2X2,i)^4z'l(rfi^2)'~^<I^^ (51) 

\ i—k ) I \kij — k i—k n—k—i 

= ^. _ (2^2,22^1,1 - 2^1,22:2,1) "2^12^1,2 

j-k 

r,k /'n-2k\ / _ ™ ^ \kS^/n-k-i\r i-k \ 

= 2 (i_J(x2,22;i,l-Xi,2X2,l) 2^( , ){j^k-s)- 

In this sum, only monomials with D = D\ ■ for some t, occur and they have coefficient 
{Anek,, e,,) 2^ ("rf ) g (-^^) (^.-^ J (^.^ J (-l)^"*-. (52) 

Note that for i ^ j the vectors ek,i and Ckj are orthogonal as they have disjoint support. 
Taking t = i = j we see that 

{ek,^,ek,^) = 2>'{\-f). (53) 
It follows that the block diagonalization is given by: 



fc=0 * '■' 
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where 



j-k 



nt I -\\j—t~s(n—2k\ \ ^ (n—k—i\( i—k \l k 

Pi,j,k ■= (-ir [i-k) s ) [j-k~s) [j-t- 



(55) 



s=0 



Remark 2. This is the same block diagonalization as was given by Schrijverf^, except that 
there a different expression (but the same value!) was used for P^j/^, namely 



E/ 1 \u—t tu\ (n—2k\ ln—k—u\ (n—k—u\ 
y-^) \t)\u-k)\ i-u )\ j-u )■ 

u=0 



(56) 



More algebras for binary codes 

Fix a positive integer t and let X := {0,1}*. The symmetric group of two elements, 5*2, 
acts on X by exchanging the symbols and 1 (in all t positions). This induces an action 
of S2 on the set oi X x X matrices by simultaneous permutation of the rows and columns. 
Let 13 := {A € C"'''^'''' | A is S'2-invariant} be the algebra of matrices invariant under this 
action. The algebra An ■= Sym"(i?) can be indentified with the set of matrices with rows 
and columns indexed by alH x n binary matrices / that are invariant under all permutations 
of the indices induced by either permuting the n rows of / or by action of S2 on any subset 
of the columns of /. 

'A 



^(C{o,i}*-ix{o,i}*-i sending 



A + B(BA-B, 



Since B ^ C^o^y-'><mV-' ^ ^ , ^ ^ 

Theorem... gives an explicit block diagonalization of An- The simplest case, t = 1 gives a 
diagonalization of the Bose-Mesner algebra of the Hamming cube. In the case t = 2 we see 
that An = ef=oSym'^(C^^^,i)Sym"~^(C^^^), a direct sum of tensor products of Terwilliger 
algebras for the Hamming scheme. 



Terwilliger algebra of the nonbinary Hamming scheme 

Let g > 3 be an integer and let B C C^^"^ be the set of matrices with rows and columns 
indexed by {0, 1, . . . ,q — 1} that are invariant under simultaneous permutation of the rows 
and columns in {1, . . . ,q — 1}. The algebra B is easily seen to have dimension 5, where 
Bi, . . . , Br, form a basis by defining 



I X3 



^XiBi :-- 



i=l 



X3 X4 X5 



X2 \ 
x^ \ 



X-i X5 



■ 2:5 

X5 X4 , 



A block diagonalization B 
(^{xxBx - 

for xi, . . . , X5 G C. 



• • • + xs^s) = 



x\ X2^q-l 



X4-XS 



(57) 



(58) 
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The algebra Aq^n ■= Sym"(S) is known as the Terwilhger algebra of the nonbinary Ham- 
ming scheme. It can be used for deriving bounds on the size of nonbinary codes from semidef- 
inite programming, see [Ij. 

Applying Theorem [HI we obtain a block diagonalization of Aq^n given by 

n 

n / , X \c n—w—a~b—c J 

\c n—w—a—b~cj ^ ' 

VD = xtx\xl{x^ + {q-2)x^r-^-^-'-\x^-x^r{q-l)"-^[x'']. (60) 
Restricting to D with yij we obtain: 



Using the block diagonalization of Sym" ""(C^^^), we obtain a block diagonalization of 
Aq^n- This agrees with the block diagonalization found in [l]. 
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